MEG-recordings of DC fields using a 151-channel whole-head device 
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1 Introduction 

Several physiological and pathophysiological phe¬ 
nomena are known to cause slow shifts of the 
electrical brain potential. Previous studies could 
also demonstrate slowly time-varying neuromagnetic 
fields (DC fields, / < 0.1 Hz) in several conditions 
(e.g., hyperventilation [1, 2], anoxic depolariza¬ 
tion [3], focal ischemia [4]). The non-invasive 
magnetoencephalography (MEG) measures electro- 
physiological parameters remotely. Therefore, it is 
not prone to impedance changes at the electrode¬ 
tissue interface like the electroencephalography, 
which hinder long-term recordings of DC shifts. 

The increase of magnetic noise in the low frequency 
domain (1// noise) interferes with the measurement 
of the neuronal DC fields. Nevertheless, recent stud¬ 
ies have shown the feasibility of DC recordings based 
on the periodic movement of the source even in a 
noisy clinical environment [2, 5]. By modulating the 
source position, i.e., the head, with respect to the de¬ 
tector, the DC magnetic field is shifted to a higher fre¬ 
quency where less environmental noise is present. 
Based on this concept, a new technique to meas¬ 
ure DC fields using whole-head magnetoencephalo¬ 
graphy (MEG) was developed. A device for moving 
the head and a continuous head localization were im¬ 
plemented and validated. In order to explore the cap¬ 
abilities of the method, artificial fields of a magnetic 
dipole phantom and biomagnetic fields were investig¬ 
ated. 

2 Methods 

All recordings were performed in a magnetically 
shielded room (Vakuumschmelze AK3b) with a 151- 
channel SQUID system (CTF Systems Inc.), using 
axial first-order gradiometers with 5 cm baseline. Its 
sensor array is shaped like a helmet, in order to min¬ 
imize the distance between the sensors and the human 
scalp and to provide whole-head coverage (Figure 1). 
The magnetic field signals were sampled at a rate of 


625 Hz. 

2.1 Mechanical modulation 

To transpose the DC fields to a higher frequency, a 
device to move the head is essential. Therefore, two 
air cushions are placed in the residual space between 
each side of the head and the detector-helmet (Fig¬ 
ure 1). The head is moved by inflating and deflat¬ 
ing the cushions alternately. The cushions are connec¬ 
ted via two flexible plastic tubes to a pneumatic sys¬ 
tem, which is located outside of the shielded room. 
The air flow is controlled by electric valves. This 
way, various modulation schemes, e.g., different amp¬ 
litudes and frequencies, are possible. In order to avoid 
additional artifacts, only non-magnetic materials are 
used inside the shielded room. 



Figure 1: Mechanical modulation: Two air cush¬ 
ions are inflated and deflated alternately to induce the 
movement of the head. Note the small gap between the 
head and detector-helmet which restricts the possible 
movements. 


2.2 Continuous head localization 

In order to extract the DC field from the modulated 
magnetic field by demodulation, the trajectory of the 



head movement must be known. Based on the head 
localization system shipped with the MEG system, we 
developed a system that monitors the head position 
continuously. Three coils are fixed to fiduciary points 
on the head (i.e., nasion, two pre-auricular points) and 
are driven with AC currents of different frequencies. 
The regression coefficients of the current of one coil 
and the sensor values reflect the magnetic field dis¬ 
tribution caused by this particular coil. Performing a 
single magnetic dipole fit on the field distribution, the 
coil is localized. From the three head coil positions 
the origin and orientation of the head is determined. 
Choosing high frequencies (> 125 Hz) for the coil cur¬ 
rents, their contribution to the signal can be removed 
by low-pass filtering before the reconstruction of the 
DC field. Hence, additional artifacts arising from the 
continuous head localization can be avoided. 

2.3 Reconstruction of DC field and sources 

To reconstruct the DC field, the recorded oscillating 
field values are demodulated by calculating the correl¬ 
ation between the field values and the monitored head 
position. Since the air cushions induce movements 
along all dimensions and even rotations, the correla¬ 
tion has to be calculated for all spatial dimensions. For 
small head movements (r max < 1 cm) and rotation 
angles (a < 2°) the DC field and its change due to the 
modulation can be described by its first-order expan¬ 
sion. Thereby, it is possible to extract the gradient of 
the DC field value at the /, -th sensor from the accord¬ 
ing correlation amplitude 

Ak= f T n(t)B k (t)dt (l) 

Jo 

exploiting 

A k = £ (n(t)rf (l)) dt VB DC ,k , (2) 

where r k denotes the relative position of the head to 
the sensor, B k the oscillating field value, and V B nc,k 
the spatial gradient of the DC field value. Obviously, 
the correlation amplitude equals the gradient of the 
DC field that is modified by the correlation matrix of 
the movement. Finally, three spatial correlation amp¬ 
litudes for each sensor and modulation period T are 
obtained. The DC field values may be reconstructed 
by numerical integration. 

To determine the sources of the DC field, a multiple 
dipole fit for both magnetic and current dipoles has 
been developed and implemented. The fit uses the 


3D spatial information of the correlation amplitudes. 
Based on a modified gradient method, it minimizes the 
residual variance of the calculated dipole model and 
measured correlation amplitudes. After determining 
the parameters of the sources, the DC field distribu¬ 
tion is estimated by forward-calculation. 

3 Results 

In order to show the capabilities of the newly de¬ 
veloped technique, both artifical and physiological 
DC fields were measured. In preliminary measure¬ 
ments the accuracy and precision of the continuous 
head localization was tested. 

3.1 Accuracy of the head localization 

The three head coils were affixed to a spherical plastic 
phantom (radius 7.5 cm) for continuos localization. 
The phantom was kept in a fixed position inside the 
detector-helmet during the measurements. Six meas¬ 
urements (duration 15 s each) with different head coil 
currents in the range of 0.05-4 //A were performed. 
Offline analysis was done with different rates of loc¬ 
alization. Since the phantom was not moved, any vari¬ 
ances in the determined positions came from localiza¬ 
tion errors. For each head coil current and localization 
rate the standard deviation was calculated. 

A nearly linear dependence of the standard deviation 
on the rate of localization within an interval of 5 Hz 
to 50 Hz was shown. The dependence on the head 
coil current was more complicated, but in general the 
precision increased with higher currents. Because a 
high localization rate and precision should be com¬ 
bined with small head coil artifacts, a medium local¬ 
ization rate of 10 Hz and a coil current in the order of 
0.4 //A was chosen for the following experiments. For 
these settings the standard deviation for determining 
the head origin was about 0.1 mm and the standard de¬ 
viation for determining the head orientation was about 
0.3°. 

3.2 Artificial DC field of a magnetic dipole 
phantom 

Three head coils were affixed to a spherical plastic 
phantom (radius 7.5 cm) for continuos localization. 
An additional coil was fixed on the surface of the 
phantom, which was driven with a DC current (I pc ~ 
2 fiA), approximating a magnetic dipole with a mag¬ 
nitude of 11m11 « 2 • 10 _9 A • m 2 . The DC field was 




Figure 2: Topography of the DC field gradient of a 
radially oriented, centrally located magnetic dipole. 
Only the component along the main movement axis 
(from left to right) is ploted. A measurement without 
the source was subtracted. The artifacts on each side 
were caused by the head coils and air cushions. 


measured for different modulation frequencies (0.5 
Hz and 1 Hz), durations (15 s and 100 s) and dipole 
positions ("central", "left temporal", and "occipital"). 
The phantom was moved harmonically with an amp¬ 
litude of approximately 5 mm along the axis connect¬ 
ing the left and right head coil. Each measurement 
was repeated with no DC current, in order to investig¬ 
ate and to be able to subtract the magnetic field of the 
modulation device and head coils. To detect possible 
systematic errors, the DC measurement was compared 
with an AC field measurement in the case of the cent¬ 
rally located dipole (Figure 2). 

After demodulation a single magnetic dipole fit was 
performed. The dipole was localized with less than 
8 mm standard deviation for all measurements. The 
localization of the AC and DC dipole differed by ap¬ 
proximately 1.2 cm, the dipole moment by less than 
5%, resembling the expected value. DC field distri¬ 
butions with a maximum peak of less than 4 pT were 
extracted from the data. Figure 2 and 3 show a topo¬ 
graphical distribution of one spatial component of the 
gradient and the DC field for the centrally located di¬ 
pole. Typical artifacts caused by the head coils and air 
cushions can be observed in Figure 2. Using the meas¬ 
urements without any DC current, the statistical error 
of the technique was estimated to be in the order of 10 
fT/m for the gradient. 



Figure 3: DC field of the same dipole as in Figure 2. 
After a single magnetic dipole fit, forward-calculation 
of the DC field was performed. 


3.3 Physiological DC field 

Preliminary measurements of neuromagnetic fields of 
two subjects (male: 26 years, female: 34 years) were 
performed, in order to show whether physiological 
DC fields could be detected. During the measurement 
the head was moved along the axis connecting both 
ears. The modulation frequency was set to 1 Hz to 
avoid discomfort. White noise was delivered via a 
plastic tube to the subject’s ear, in order to suppress 
evoked potentials caused by the ventilation noise of 
the air cushions. The subject fixated a target to reduce 
eye movements. To monitor the head position, three 
head coils were attached to the head at three fiduciary 
points (i.e., nasion, pre-auricular points). 

The subject’s head was moved nearly harmonically 
with an amplitude of about 4 mm. The rotations of the 
head were below 2°. 

To elicit neuronal DC fields, the subject put her fin¬ 
ger (digit IV) on a periodically moving nylon brush. 
Alternately, the left, right or no finger was stimulated 
for 20 s. This sequence was repeated in different or¬ 
der three times. Altogether six experiments were per¬ 
formed. 

Figure 4 shows an example of the gradient of the DC 
field. In the frontal area large gradients occured, pos¬ 
sibly caused by the physiological DC field of the eyes. 
The magnitude and topographical pattern of the gradi¬ 
ent is in agreement with models, which describe the 
magneto-oculagram (MOG) produced by eye move¬ 
ments and blinks [6, 7]. 

The source of the residual noise is not known exactly. 
We suspect contamination of the hair, skin and equip¬ 
ment with magnetic substances (e.g., artifacts arising 
from adhesive tapes were observed). 
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Figure 4: Typical topography of the DC field gradient 
reconstructed from the physiological data. Only the 
component along the main movement axis (from left to 
right) is ploted. The large field gradient in the frontal 
area was most likely caused by the DC field of the sub¬ 
ject’s eyes. 


No task-related neuronal DC field could be extracted 
from the data. Spontaneous neuromagnetic fields and 
the artifacts might have masked possible task-related 
DC fields. The physiological effect could have been 
too small, as well. 

4 Discussion 

A new technique to measure DC fields using whole- 
head magneotencephalography was developed. 
Based on mechanical modulation with air cushions 
the movement of the human head inside the detector- 
helmet could be induced successfully (amplitude « 4 
mm). Continuous monitoring of the movement was 
possible using AC current driven head coils. An ap¬ 
propriate localization rate of 10 Hz and a localization 
accuracy of 0.1 mm and 0.3° were achieved. Both 
physiological and artificial DC fields were measured. 
DC fields with a maximum peak of less than 4 pT 
were extracted. 

Yet, different sources of artifacts and field distortions 
have to be considered and further increase of the sens¬ 
itivity is needed to detect neuromagnetic fields. Cur¬ 
rently, we are developing methods to degauss mag¬ 
netic contaminations, which are limited by safety 
rules for the exposure to electromagnetic fields [8]. 
The modulation technique using whole-head magne¬ 
toencephalography offers new perspectives to mon¬ 
itor intracranial DC currents (e.g., lateralization ef¬ 
fects). 
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